Modeling Arcs 



Z. Insepov, J. Norem* 
Argonne National Laboratory, Argonne, IL 60439, USA 

S. Veitzer, S. Mahalingam 
Tech-X, Corp., Boulder CO 8080^ 

Although vacuum arcs were first identified over 110 years ago, they are not yet well understood. 
We have since developed a model of breakdown and gradient limits that tries to explain, in a 
self-consistent way: arc triggering, plasma initiation, plasma evolution, surface damage and gra- 
dient limits. We use simple PIC codes for modeling plasmas, molecular dynamics for modeling 
surface breakdown, and surface damage, and mesoscale surface thermodynamics and finite element 
electrostatic codes for to evaluate surface properties. Since any given experiment seems to have 
more variables than data points, we have tried to consider a wide variety of arcing (rf structures, 
e beam welding, laser ablation, etc.) to help constrain the problem, and concentrate on common 
mechanisms. While the mechanisms can be comparatively simple, modeling can be challenging. 



I. INTRODUCTION 

The field of vacuum arcs is very old [T-^. It is complicated by the large number of variables, their large parameter 
range and the number of mechanisms involved, compared to the limited range of measurements that are done on 
individual arcs. This work began with the goal of understanding the production of dark currents in rf cavities, where 
X ray production would interfere with a planned particle beam experiment [IHH]. We found, however, that by cleanly 
looking at the field emission from asperities at the breakdown limit, we could directly study the environment of the 
pre-breakdown surface. Initial effort produced a model of the breakdown trigger based on Coulomb explosions, and 
subsequent work has examined the properties of arcs, primarily by assuming that unipolar arcs are the dominant 
mechanism [9HTT]. 

II. THE ARC MODEL 

We divide the problem of vacuum arcs into four parts: the trigger, plasma initiation, the plasma growth phase, and 
plasma damage of the metallic surface. At all stages the driving mechanism seems to be the plasma and surface 
electric fields, and Ohmic heating of the surface is not required either by the trigger or for plasma evolution. We 
believe that Coulomb explosions and unipolar arcs can explain much of the behavior of these arcs. 

A. Triggers 

Our work measuring dark currents and x rays from rf cavities showed that rf structures operate in a mode where 
the dark currents / depend on the surface electric field E like / ~ E^^ — E'^^ independent of state of conditioning, 
frequency or stored energy. When the Fowler-Nordheim currents are plotted in this way we find that this behavior 
is characteristic of field emission from a surface at a field of ~10 GV/m. At this field the electric tensile stress. 
(J = eoE'^/2 (eo being the permittivity of free space) from the electric field is comparable to the tensile strength of 
the material, and the material is subject to fatigue failure due to the rf oscillations. 

While the conventional wisdom is that arc triggers are due to Ohmic heating in the surface, asperities of the required 
dimensions have not been found, and arc triggers occur randomly, with no time to warm up. Coulomb explosions, 
augmented by fatigue, do not significantly constrain the breakdown site geometry. Note that while data shows that 
the breakdown rate is proportional to ~ E^^^ this behavior is compatible with a) Ohmic heating, b) failure due fatigue 
stress and c) electromigration, so is inconclusive. 
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FIG. 1: a) Measurements of dark currents from a 1.3 GHz accelerator cavity showing comparisons with Fowler- Nor dheim 
theory, with enhancement factors and emitter area factors 4 . b) Determining the surface field from the exponent / ~ E^, for 
a work function of 4 eV. c) Molecular Dynamics simulation of Maxwell stresses pulling copper asperities apart, which would 
trigger the arc 6 . 



B. Plasma Initiation 



The defining property of a vacuum arc is that it can occur in a vacuum. Many experiments over the years have 
shown that the arc can occur at very large gaps, implying a single surface phenomenon [3]. These conditions imply 
that field emitted beams must ionize material fractured off the surface by electric tensile stress. In order to efficiently 
ionize material close to the cathode, our OOPIC PRO simulations show that this requires a total mass of material 
near the field emitter equivalent to half of a monolayer [12 . This density of gas is sufficient to accommodate a plasma 
density of 10^^ m~^. 

The ions produced are driven from their production point by space charge fields, that produce a streaming, almost 
mono-energetic, ion current that hits the surface producing a variety of effects. 



C. Plasma Evolution 



Once some plasma exist above the surface the surface electric field will be affected by both the applied electric field 
and the surface field produced by the sheath [13 . As the density, n, increases, simulations show that a number of 
things happen: 1) the sheath potential, (j) stays roughly constant, since this potential is a function of the overall net 
charge, which evidently remains constant, 2) the surface electric field increases diS E ^ (jy/Xo ~ 0V^5 where Xd '^s the 
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FIG. 2: Time evolution of rf breakdown at 805 MHz and the temperature (left) , and density of the plasma as calculated by 
OOPIC Pro (right). 

Debye length, 3) the increase in the surface electric field further increases the field emission current and the number 
of emitters. We find that the codes show exponentially increasing arc densities, and believe that the arc eventually 
becomes a non-Debye plasma. Our model the development of rf arcs is shown in Fig 2. We expect that DC arcs 
would follow a similar path in parameter space. 

OOPIC Pro simulations of show that the sheath potential is on the order of 75 volts during field emission, however 
this is evaluated for a small, roughly spherical plasma in contact with the surface and is due to the local ionization. 
This potential decreases in all directions, however, so the model should be consistent with lower burn voltages seen 
in small gap arcs. 

In accelerators there seem to be two limiting cases for arc interactions with their environment, which we have 
described as killer arcs and parasitic arcs, where the killer arcs are able to directly short out the driving potential and 
extinguish themselves, and parasitic arcs are unable to dispose of sufficient energy to perturb the driving field, and 
can survive for as long as external potentials are available. 

We find that the growth time of the plasma are consistent with measurements we can make of fully formed arcs 
using X rays, however, the computer simulations stop before experimentally measurable currents are produced, so 
some extrapolation is required. 

D. Surface Damage 

Surface damage by arcs can take many forms, depending on the parameters of the plasma and the properties of 
the surface. We find that the dominant interaction is self sputtering by plasma ions streaming out of the plasma. As 
the plasma evolves it becomes denser and the Debye length shortens, increasing the surface field and field emission 
currents and ion flux (Fig lb.). We have calculated the self sputtering rates for high surface fields, high surface 
temperatures and varying grain orientation with the results shown in Fig 2. 

Although the temperature of the plasma is very low, the potential induced by ionization produces a very high plasma 
pressure for the high density plasma. After the surface is melted, the plasma pressure p = nkT, or nkE^ pushes on 
the surface, the electrostatic field pulls on the plasma and the surface tension tends to flatten the surface. The 
combination of a pulling (or pushing ) force on the surface, combined with the surface tension produces an instability 
in an initially fiat surface which results in a spinodal decomposition - the formation of ripples in the surface. The 
dimensions of these ripples can then be used to estimate the parameters of the plasma. Ripples can also be produced 
by an oblique flow of ions carrying momentum to the surface in the same say as wind can produce ocean waves. Thus 
ripples of two types can be produced on melted surface. 

The high plasma pressure can produce sufficient liquid motion to generate particulates with sufficient mass and 
volume to uniformly cover the surface with secondary breakdown sites. 



4 




Temperature, K 



FIG. 3: Self sustained arcs require high self sputtering yields. Calculations show either high surface fields (left) or surface 
temperatures (right) will provide these high yields. 



1 . Self- sputtering 

The intense flux of low energy ions hitting the surface is the primary mechanism by which the plasma affects the 
surface. In addition the deposited heat and pressure, self-sputtering helps to determine the evolution of the plasma 
by determining the atomic fluxes back to the plasma. Anders has shown that self-sputtering coefficients significantly 
above 1 (we assume 10) can produce a self sustaining arc [11 . Since the surface environment beneath the plasma 
is poorly understood, we have used Molecular Dynamics to estimate the self-sputtering coefficients, as a function of 
surface temperature, surface electric field and grain orientation. 

• The effect of temperature is the simplest. As the temperature approaches the melting point, the surface binding 
energy decreases and the sputtering yield increases. Because surface atoms are more mobile slightly below the 
melting point, this increase occurs slightly below the bulk melting temperature as shown in Fig 3b. 

• The tensile stress induced by a strong electric field tends to pull atoms out of a surface. In field evaporation 
at fields of ~30 GV/m, this field is sufficient to pull individual atoms out of a polished surface. We find, using 
Molecular Dynamics, that fields on the order of 1 - 3 GV/m, which we expect beneath plasma, are sufficient to 
significantly modify the low temperature self sputtering coefficient. 

• Surface grain orientation also seems to also affect self-sputtering yields. Recent studies of rf systems has shown 
surface damage that seems to be dependent on grain orientation. We have also calculated the sensitivity of the 
self-sputtering coefficient on grain orientation. We find that this effect can be quite significant at low energies. 



2. field enhancements 

Although we, like others, do not find whiskers on our surfaces, we do find many cracked regions and many sharp 
edged craters. Fig 3 shows an example of one array of cracks consistent with the surface cooling by about 1000 deg. 
When we model the field pattern on crack junctions, we find that the crack junctions have a field enhancement on 
the order of 100, which is consistent with the values we measured from dark currents. This enhancement would be 
multiplied by a factor determined by the local shape of the surface, so field enhancement as large as 1000 are in 
principle possible. In this model the area of the emitters would be only a few nm^, but a large number of them 
should be produced. This analysis is described in Fig 3. Data showing large area emitters can be accommodated by 
assuming that many emitters are detected. Note that if the current is proportional to E^^^ small systematic errors in 
measuring the electric field would produce much larger errors in measuring the emitter area, particularly if a proper 
two dimensional (/3, area) fit of the data was made. 

The conventional wisdom is that field emitters heat, vaporize and ionize due to Ohmic heating due to field emission 
currents. With corners or conical geometries, the field emitting region is very small (nm^) and the thermal diffusion 
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FIG. 4: a) SEM image of sub micron cracks, b) Computer model of these crack junctions, c) Enhancement factors of electric 
field at the tip of these crack corners. 




EnhMMnKni Fa«^ 

FIG. 5: Numerical simulation of heating of a blunt corner during field emission, a). Because of the small heated volume and 
large thermal diffusion volume, there is little temperature fluctuation. The spectrum of enhancement factors from a number of 
different sources, b) . The spectrum is seems to show a e~^^^ ^ dependence in many data sets 

volume for nanosecond timescales is much larger (/im^), so the heat is quickly removed and there is no significant 
temperature increase. This effectively prohibits Ohmic heating as a breakdown trigger in many cases. 

III. GRADIENT LIMITS 

We find that the spectrum of field enhancements for many different kinds of experiments has an exponential form 
with density of sites n(/3) ~ e~^l^ . If we assume that the total number and spectrum of particulates (breakdown 
sites) is proportional to the energy in the arc, it is possible to argue that and equilibrium enhancement factor for a 
given system will be produced where the equilibrium enhancement factor /3eg ~ In [/, where U is the energy of the 
arc. We have shown that this behavior is consistent with scaling laws such as the Kilpatrick limit. 

Arcing phenomena are a special case of plasma surface interactions a field under active study for many reasons. We 
believe that breakdown in high gradient accelerators, surface defects from e beam welding, laser ablation, tokamak rf 
limits, micrometeorite impacts, and some some power grid failure modes should have many common mechanisms and 
is should be useful to identify and study them more systematically. 
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IV. COMMENTS 

The study of arcs is complicated because the number of mechanisms involved in any experiment usually exceeds the 
number of experimental variables. We believe that the basic assumptions of arc models should be carefully examined 
as the conventional wisdom over-constrains the problem. For example, starting with Dyke, et al., the conventional 
wisdom is that breakdown is initiated from ohmic heating of whiskers or asperities of comparable geometries. The 
whisker geometry is unique, because the ratio of volumes in which Ohmic heating and thermal diffusion occur are 
comparable. With all other geometries heating must be much less, and a different mechanism of plasma production 
like Coulomb explosions, more compatible with the shapes that are actually seen in SEM images, seems to be required. 

A. R&D issues 

There are a number of general questions. Although arcing occurs in many different fields, it is interesting to see 
how much of arc behavior is common to all arcing applications. Beyond that, there are a number of more specialized 
questions that may become accessible. How do non-Debye plasmas interact with materials? What are the damage 
mechanisms? How do the presence of preexisting plasmas and strong external magnetic fields in three dimensions 
affect the triggering and evolution of an arc? How are the equilibrium field enhancements, and gradient limit affected 
by these variables? What is the effect of dense neutral gas on the arc? To what extent can surface treatment mitigate 
or modify gradient limits? What determines the time structure of these arcs? These questions will test modeling. 

V. CONCLUSIONS 

We have developed and modeled a picture of rf arcs in accelerators, where the basic mechanisms are Coulomb 
explosions and unipolar arcs. One of the problems with this field, however, is that there frequently are alternative 
mechanisms that can produce satisfactory comparisons to a given set of data. We believe it is important to both 
increase the precision and self-consistency of the models, but it is important to extend models to the widest range of 
phenomena as the most useful test of their value. 
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